Two-dimensional electron liquid state at LaA103-SrTi03 interfaces 



O 

(N 



c/5 



c 

o 

> 



o 
o 



X 



M. Breitschaft, 1 V. Tinkl, 1 N. Pavlenko, 1 - 2 S. Paetel, 1 C. Richter, 1 J. R. 

Kirtley, 1 Y. C. Liao, 1 G. Hammerl, 1 V. Eyert, 1 T. Kopp, 1 and J. Mannhart^E 

Center for Electronic Correlations and Magnetism, 
University of Augsburg, D-86135 Augsburg, Germany 
2 Institute for Condensed Matter Physics, 1 Svientsitsky Str., UA-79011 Lviv, Ukraine 

(Dated: May 5, 2010) 

Using tunneling spectroscopy we have measured the spectral density of states of the mobile, two- 
dimensional electron system generated at the LaAlC>3-SrTi03 interface. As shown by the density 
of states the interface electron system differs qualitatively, first, from the electron systems of the 
materials defining the interface and, second, from the two-dimensional electron gases formed at 
interfaces between conventional semiconductors. 



Two-dimensional (2D) conducting electron systems are 
generated at interfaces between a large variety of insu- 
lating oxides^. These interfaces show a broad spec- 
trum of different properties. The quantum-Hall effect 
has been found 3 -, for example, for the electron system at 
the ZnO-(Mg a; Zni_ a: 0) interface. For the interface be- 
tween LaTiC>3 and SrTiOs, first explored experimentally 
by Ohtomo and Hwang 4 , spin ordering and ferro-orbital 
ordering has been predicted^. The most widely investi- 
gated electron system at oxide interfaces is the metallic 
state created at the interface between the charge transfer- 
insulators LaA103 and Ti02-terminated SrTiO^ . This 
electron system forms a 2D superconductor with a T c of 
R3 250 mK that is easily tunable by electric gate fields 6 . 
The LaA103-SrTiC>3 interface has also been predicted^ 
and reported 8 - to develop magnetic order. 

While several theoretical models have been devel- 
oped to describe the electronic properties of these 
interfaces 9 - - — , less information on the electronic struc- 
ture has been provided by experiments. At room tem- 
perature its thickness has been inferred from scanning 
tunneling microscopy (STM) writing experiments^ 3 -, from 
photocmission 14 , and from cross- sectional STM (Ref.[l]| 
to be at most a few nanometers. Further, hard x-ray pho- 
toelectron emission has shown that the charge carriers 
at the interface occupy Ti3d states^. Studies of x-ray 
absorption spectroscopy furthermore revealed that ener- 
getically the crystal-field split Ti levels are rearranged, 
such that the id xy levels are the first available states for 
the conducting electrons 1 ^. 

The spectral density of states (DOS) at the interface 
is a fundamental property that characterizes the elec- 
tron system. As it furthermore can be calculated as 
well as measured, it is a key property for the under- 
standing of the electron system at the interface. For 
measurements of the spectral DOS scanning tunneling 
spectroscopy (STS) is a powerful technique^, which has 
been used extensively to characterize 2D electron gases 
(2DEGs) in semiconductor systems. STS was employed, 
in particular, to probe surfaces of semiconducting thin 
films where electrons arc confined by the film thickness 18 . 
STS was also used successfully to analyze cross-sectional 
cleavage planes of semiconductor heterostructure o 19 ' 20 . 



In addition, semiconductor surfaces, below which elec- 
trons are confined in band bending regions induced by 
ion implantation^ or surfaces at which electron gases 
were generated by adsorbates22, were explored. 

Here we report on STS measurements of the spectral 
DOS of the electron system at the LaA103-SrTi03 in- 
terface. We find the measured DOS to be in excellent 
agreement with the interface DOS calculated in density- 
functional theory (DFT), providing evidence that the 
tunneling current in the STS measurements is carried by 
interface states. The measured spectrum of the interface 
DOS and therefore the electron system differs qualita- 
tively from the DOS of doped bulk SrTi0 3 or LaA10 3 . 
The electron system cannot be accurately described as 
a thin layer of doped SrTiOs. The measurements reveal 
furthermore that the electron system also differs quali- 
tatively from the hitherto known 2D electron systems at 
interfaces between conventional semiconductors. We find 
the electrons confined in multiple layers of quantum wells 
given by the ionic potentials of the Ti06 octahedra. In 
these wells the electrons are subject to the correlations 
characteristic of the d orbitals of the Ti ions. The spec- 
tral DOS is not a step function as is the case for standard 
semiconductor interfaces but rather resembles the DOS 
of Ti3d states. Quantum wells and electronic systems of 
this kind are unknown from the 2DEGs in conventional 
semiconductors, in graphene or in ZnO. 

For the studies, we fabricated LaA103-SrTi03 het- 
erostructures with 4 unit-cell (uc) thick (« 1.6 nm) epi- 
taxial LaA103 layers to obtain measurable tunneling 
currents. This thickness was chosen because it is the 
minimum thickness required to generate the conducting 
interface 23 . For larger LaA10 3 thicknesses the tunnel- 
ing current densities become impractically small. The 
samples were grown by standard pulsed laser deposition 
as described in Ref. [2J. For deposition the SrTiOs sub- 
strates were heated to 780 °C in an oxygen background 
pressure of 8 x 10~ 5 mbar. The LaA103 film growth was 
monitored by reflection high-energy electron diffraction. 
While SrTiOs surfaces are known to show numerous sur- 
face reconstructions^—, x-ray diffraction showed no ev- 
idence of distortions of the LaA10 3 films, which could 
be attributed to a SrTi0 3 surface reconstruction, sug- 



2 




4 unit cells LaAI0 3 
interface electron system 
SrTi0 3 



FIG. 1: Illustration of the experimental configuration. A 
metallic state is formed at the interface between the SrTi03 
substrate and a 4 uc thick layer of LaA103. Scanning tunnel- 
ing microscopy and spectroscopy are performed by monitoring 
the tunneling current It between the tip and the sample as a 
function of V s , the voltage at the sample relative to the tip. 
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gesting that the LaAlC>3 growth stabilizes the standard 
SrTiOs structure at the interface. Titanium plugs filling 
ion etched holes were used to contact the interfaces. After 
a heating procedure in a preparation chamber—, the sam- 
ples were transferred in situ to the scanning probe mi- 
croscope (SPM), which operates in ultrahigh vacuum at 
4.7 K. An iridium spall attached to a cantilever based on a 
quartz tuning fork 33 with a spring constant of 1800 N/m 
was used as a tip. The tip was treated in situ by field 
emission 24 . The cantilever was not excited mechanically 
during STM and STS measurements. The experimental 
setup is sketched in Fig. Q] Typical measurement pa- 
rameters were tunneling currents of 10 pA, sweep rates 
of 0.01 V/s and scanning speeds of lOnm/s. 

Imaging the LaA103-SrTiC>3 heterostructures by fre- 
quency modulation scanning force microscopy 3 ^ (FM- 
SFM) as well as by constant current STM revealed the 
standard step-and-terrace structure resulting from the 
slight miscut of the SrTi0 3 substrates (Fig. [2]). While 
on more conventional samples excellent resolution was 
achieved with the SPM employed 3 ^, it was impossi- 
ble to obtain atomic resolution on the LaA103-SrTi03 
heterostructures^. 

Conductance-voltage characteristics dI t /dV s (V s ) were 
measured using a standard lock-in technique 2 ^. Simulta- 
neously, the tunneling current was measured as a func- 
tion of voltage. The normalized differential conductance 
NDC = (dI t /dV s )/(It/V s + e) was determined as a mea- 
sure of the sample DOS 3 —. The spectra were taken on 
sample areas where the step-and-terrace structure was 
resolved in STM topography. We found the characteris- 
tic spectroscopic features to be reproducible across four 
samples 24 . Figure G^a) shows a representative depen- 
dence of the NDC on voltage. The conductances are 
minute for negative voltages (tunneling from occupied 
sample states). For positive voltages (tunneling into un- 
occupied sample states) the spectroscopically accessible 



FIG. 2: Scanning probe microscopy images of LaA103-SrTi03 
heterostructures. (a) Topographic FM-SFM image of the 
LaA103 film, (b) Profile taken along the line indicated in 
(a), (c) Topographic STM image acquired on the LaA103- 
SrTiOs heterostructure recorded with a scanning speed of 
5nm/s, a bias voltage V s = 2 V, and a current set point of 
10 pA. (d) Profile taken along the line indicated in (c). All 
data were taken at 4.7 K. For further detail see Ref. |24|. 

energy range is limited at low voltages by small tunneling 
conductances and at high voltages by large electric fields 
destabilizing the tunneling gap. To measure the tunnel- 
ing characteristics at a given sample location over a large 
voltage range, several spectra were therefore taken at dif- 
ferent tip-sample separations as determined according to 
the tunneling currents it, stab at given gap voltages V^ s t a b- 
Three characteristics measured with different tip-sample 
separations, from two different samples, are shown in Fig. 
[J2a). In these spectra, clear peaks are seen at « 0.6, 
« 0.8, 1.4, and « 1.8 V. 

To identify the electron states carrying the measured 
tunneling current and to explore the role of electronic 
correlations at the interface we compare the measured 
DOS to predictions of DFT. We performed local-density 
approximation (LDA) and LDA+U calculation s 37 ' 38 of 
the layer-resolved DOS of LaA103-SrTi03 heterostruc- 
tures. Further information on these calculations is given 
in Ref. HJ. While differences are present in details, the 
calculated state densities and the effective electron mass 
of 3.25 bare electron masses are consistent with those 
reported in Refs. iHH! and 53 

In Fig. Mb) the Ti 3d DOS of the interface Ti0 2 layer 
calculated using LDA for a supercell with a 4 uc thick 
LaA103 layer on SrTiOs is shown. According to the cal- 
culation, electronic reconstruction leads to a doping of 
O 2p states located in the topmost AIO2 layer with holes 
and of Ti 3d i 2 g states located at the interface with elec- 
trons. Experimentally, it is only the interface which is 
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FIG. 4: Illustration of the configuration of 2D electron sys- 
tems in standard semiconductor interfaces and at the LaA103- 
SrTiOa interface, (a) At the interface between the semicon- 
ductors an electron gas is generated at a potential well cre- 
ated by band bending, which typically has a width of tens of 
nanometers as determined by the electronic screening length 
s. The electron states can be approximated by the states of 
free electrons in this potential well, (b) At the oxide interface 
the potential well is provided by the Coulomb potential of the 
titanium ions in the TiOe octahedra and, to a smaller extent, 
by band bending. These potential wells are narrower than 
those at semiconductor interfaces; the "resonant" electron 
states are well approximated by the Ti3d£2 g states, which 
form a 2D electron system extended parallel to the interface. 
Due to the electronic correlations of the oxide lattices, the 
mobile electrons form an electron liquid. 



FIG. 3: Comparison of the measured differential conduc- 
tance and DFT calculated state densities, (a) NDC(V4) = 
(dIt/dV B )/(It/V B + e)(K) characteristics with e = lpA/V 
measured at several sites located far away from topographic 
steps. The measurements were performed at 4.7 K with fixed 
tip-sample separations. The different colors reflect different 
tip-sample separations (purple: it.stab = 70 pA, V s , st ab = 
2.4 V; blue: 7 t , stab = 12 pA, V S!Stab = IV; red: 7 t , stab = 12 pA, 
V SjS tab = 0.8 V). The data were averaged over an interval of 
75 mV. The purple characteristic was measured on a different 
sample than the red and blue ones, (b) Ti 3d DOS of the in- 
terface Ti02 layer calculated using LDA for a supercell with 
a 4 uc thick LaAlOs layer, (c) Ti3d DOS of the interface 
Ti02 layer calculated using LDA+C/. In (a)-(c) the positions 
of characteristic features in the NDC are marked with gray 
lines. 



found to be conducting. In the total DOS between and 
2eV the Ti3<ii2g orbitals located in the interfacial TiC>2 
layer prevail. The other, small, contributions are pro- 
vided by the TiC>2 planes of adjacent SrTiC"3 layers and, 
below ?» 0.5 eV, by the 2p states of the surface. The 
Ti3de g states contribute at energies above ps 2.8 eV and 
the La5d states at energies above « 2.2 eV. The mea- 
sured peaks at « 0.8, « 1.4, and ps 1.8 V are also present 
in the calculated DOS. The measured peaks at ps 0.6 and 
~ 1 V, however, are not represented in the LDA result. 

The LDA+C/ calculations of the interface electron sys- 
tem consider an on-site Coulomb repulsion U — 2 eV and 



a Hund coupling J = 0.8 eV in the Ti 3d shell. The choice 
of conservative values for U and J does not imply these 
values characterize the system best. Figure EJc) shows 
the DOS calculated for the supercell using LDA+C/. The 
DOS exhibits additional peaks at ps 0.6 and ps 1 eV, gen- 
erated by the splitting of the Ti 3d xz + 3d yz bands due to 
the interorbital interactions caused by the finite U and 
J. Remarkably, these peaks are observed experimentally 
but are missing in the LDA DOS. 

We note that the experimental hump at ps 1.8 eV is 
broader than the corresponding structure in LDA. How- 
ever, LDA+C/ generates a structure of approximately 
the measured width but with finer structures. These 
fine structures reflect the formation of the upper Hub- 
bard bands, which is a fundamental effect of correlated 
electron systems, arising when U is on the order of the 
bandwidth or larger. Indeed, the calculated width of the 
Ti3dt 2g band is ps 2eV= U. 

The good agreement between experiment and calcula- 
tion suggests that the electron states carrying the mea- 
sured tunneling current are the ones calculated in DFT. 
For energies between 0.5 and 2 eV, mainly Ti3d xz + 3d yz 
and Ti 3d xy orbitals of the interface Ti02 layer contribute 
to the calculated DOS and the prominent peaks result 
from these orbitals; tunneling occurs into Ti 3d orbitals 
at the interface, the significant contributions arising from 
the Ti3cit2g states. These results are consistent with the 
results of recent photoabsorption measurements^, from 
which it was concluded that the lowest unoccupied states 
are Ti 3d xv states. 
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The fact that the experimental DOS is matched sig- 
nificantly better by the LDA+C/ calculation than by the 
LDA calculation provides evidence that the electron sys- 
tem at the interface is correlated with substantial values 
of U and J on the Ti 3d orbitals. Therefore, this 2D elec- 
tron system is a liquid. This electron liquid is formed 
by correlated electrons, which can move parallel to the 
interface, but are constrained in their perpendicular mo- 
tion by the Coulomb potentials of the titanium ions of 
the final TiC>2 layers and also, to a smaller degree, by 
band bending (Fig. HJ). 

Interfaces in oxides therefore broaden the spectrum of 
available 2D electron systems from the 2DEGs of con- 
ventional semiconductors to also include systems with 
sizable electronic correlations. Such correlation effects, 



in combination with the already intriguing physics of 2D 
electron systems, promise unprecedented electronic phe- 
nomena. Influenced by electronic correlations generated 
in the ionic lattices of the oxides, electron systems at 
oxide interfaces have exceptional properties, possibly en- 
abling devices with hitherto unknown characteristics. 
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A. Sample preparation 

The LaAlOs films were epitaxially grown by pulsed laser deposition on TiCVterminated 
(001) surfaces^ of SrTi0 3 single crystals. The deposition was conducted at a substrate 
temperature of 780 °C and at an oxygen background pressure of 8 x 10~ 5 mbar. The film 
growth was monitored by reflection high-energy electron diffraction. The samples were 
cooled in 400 mbar of oxygen. During the cooling the samples were held at 600 °C for one 
hour. After cool-down an aluminum shadow mask was applied to define a stripe along one 
edge of a sample. Within this stripe the SrTiC>3 substrates were laid bare by Ar ion etching 
and then covered by electron beam evaporation with titanium to provide a contact to the 
interface electron system. 

B. Sample treatment by heating 

To achieve reproducible surface conditions after the samples were transferred in air to 
the preparation chamber of the SPM, they were radiatively heated to ~ 400 °C in an oxygen 
pressure of 10~ 2 mbar (background pressure ~ 5 x 10~ 7 mbar). 

C. Tip treatment 

A tip treatment by field emission at gap voltages up to = 400 V and emission currents 
up to several microamperes, while the tip was placed in proximity to the LaA103 surface, was 
necessary to obtain reproducible results in topographic STM imaging and STS. Usually, the 
field emission procedure damaged the sample locally. We do not expect this tip preparation 
technique to yield a clean metallic tip. The remaining contamination made it difficult to 
establish stable tunneling. The tunneling current I t and the differential conductance dI t /dV s 
tended to fluctuate on timescales from fractions of a second to several minutes, which we 
attribute to spontaneous changes of the tunneling gap configuration. By repeating the tip 
preparation procedure described, however, it was possible to reduce the fluctuations such 
that reproducible experimental results were obtained. 
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D. FM-SFM imaging (Fig. 2(a) and (b)) 



The image in Fig. 2(a) was recorded with a scanning speed of 16nm/s and a bias voltage 
of V s = 0.1 V. The eigenfrequency of the cantilever was f = 24296.8 Hz, the quality factor 
Q = 17832 and the set oscillation amplitude A = 1.8 A. The step heights in Fig. 2(a) were 
used to calibrate the scanner of the SPM. As shown by Fig. 2(a) and (b), in FM-SFM the 
surfaces are flat with an apparent roughness of 0.7 A rms . This value does not characterize the 
surface morphology quantitatively, because dissipation of vibrational energy of the cantilever 
may have caused topographic artifacts. 

E. STM imaging (Fig. 2(c) and (d)) 

We attribute the apparent surface roughness of 1 A rms seen in STM topography (Fig. 
2(c) and (d)) mainly to local variations of the tunneling barrier height, likely induced by 
adsorbates. 

F. STS lock-in technique 

The differential conductance dI t /dV s was measured by adding to the bias voltage V s an 
AC voltage with an amplitude Vmod of several millivolts and a frequency / mo d on the or- 
der of 100 Hz. The resulting in-phase AC component of the current signal was measured 
using lock-in detection. The ratio of its amplitude AJ t and Vmod approximates the differ- 
ential conductance, AI t /V mo d ~ dI t /dV s . Differences between this measured differential 
conductance and numerical derivatives of current-voltage characteristics were found to be 
insignificant. Tunneling current and differential conductance were recorded simultaneously 
with an acquisition time of 120 s per sweep direction. 

G. STS parameters in Fig. 3(a) 

All spectra in Fig. 3(a) were measured with V mo d = 5mV rms and / mo d = 77 Hz. 
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H. Reproducibility of STS measurements 

Contaminants located between the tip and the LaA103 surface may possibly distort the 
spectroscopic results. To understand the role of contaminants we investigated four samples. 
We regularly obtained stable tunneling for gap voltages up to = 1.2 V and observed the 
three characteristic peaks between V s = 0.5 V and V s = 1.2 V shown by Fig. 3(a). As we must 
expect that in these studies the configurations of adsorbates in the tunneling barrier were 
not identical from experiment to experiment, we conclude that the spectroscopic features 
cannot be caused by contaminants. 

However, adsorbates may induce fluctuations frequently observed in the tunneling current 
I t and the differential conductance dI t /dV s . The spectroscopic resolution was limited in most 
cases by such fluctuations. As explained earlier, a tip treatment by field emission made it 
possible to obtain stable tunneling, yielding reproducible spectra with low noise and low 
hysteresis. 

As shown in Fig. IST] the characteristic features of the spectra (Fig. 3(a)) were repro- 
ducible across several samples. In Fig. ISTT g) and (h) characteristics of the NDC are shown, 
which were measured at the same sample site and at nominally equal tip-sample distances. 
After the acquisition of the spectrum in Fig. ISTT g). the tunneling gap configuration ap- 
parently changed. The measurements following the acquisition of the characteristic in Fig. 
ISTT g) were accompanied by an increasing destabilization of the tunneling conditions, result- 
ing in relatively large fluctuations of the tunneling current and differential conductance. The 
spectroscopic peak structure was partially obscured. During further measurements the fluc- 
tuations of the tunneling current and differential conductance decreased, which we attribute 
to the formation of a new, stable tunneling gap configuration. The characteristic shown in 
Fig. [STT li) measured with this gap configuration differs in overall shape from the characteris- 
tic in Fig. ISITf g). In both spectra, however, the characteristic spectroscopic features appear 
at almost the same voltages, independent of the tunneling gap configuration. 

For gap voltages V s larger than 1.5 V stable tunneling was rarely obtained. Characteristics 
with a high resolution between V^, = 1.5 V and V B = 2 V as shown in Fig. 3(a) and Fig. 
ISTT g). (h) were obtained during one measurement session only. 

All tunneling spectra shown in this work were acquired at moderate tip-sample distances, 
avoiding strong repulsive force interactions between the tip and the sample which could 
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possibly induce mechanical strain on the sample. Within this distance regime, tunneling 
conductances became minute at gap voltages below 0.5 V and the signal-to-noise ratio of 
the tunneling current became small. Occasionally, it was nevertheless possible to derive the 
sample DOS for energies between OeV and 0.5 eV. Figure ISITa) and (b) compare two spectra 
where this seems to be the case. Albeit relatively high hysteresis, both NDC characteristics 
share distinct features. 

I. Possible contribution of O 2p states at energies below 0.5 eV 

According to the DFT calculations, between 0.5 eV and 2eV Ti 3<i states dominate the 
DOS. For smaller energies also 02p states of the topmost AIO2 layer adjoining vacuum 
contribute and may play a role in tunneling. Figure [S2] compares the measured differential 
conductance NDC (Fig. IS2(a)) with the 2p DOS of the topmost AIO2 layer and the Ti3d 
DOS of the interface Ti02 layer calculated using LDA+U (Fig. IS27b)). 

J. LDA and LDA+C/ calculations 

The spectral density of states of the LaA103-SrTi03 bilayer was calculated for a supercell 
comprising 4 uc LaA103, SrTiOs layers (~ lnm), 4 uc LaAlOs and vacuum (^ 1.3 nm), as 
sketched in Fig. [S3] We employ the WIEN2k program package^, using 120 fc-points (15 
fc-points in the irreducible wedge of the first Brillouin zone) for the LDA calculation, 600 k- 
points (21 fc-points in the irreducible wedge of the first Brillouin zone), a Coulomb repulsion 
U = 2eV and a Hund coupling J = 0.8 eV in the Ti3e? shell for the LDA+C/ calculation^. 
To avoid a spurious mixing of the La / states with the Ti3d bands, a large U of 8eV is 
imposed on the La / states, a procedure which was introduced in Ref. 0. The calculations 
involved structural relaxation along all crystal axes. 

Figure IS4l shows the total DOS of the full supercell (green), the DOS of the interface Ti02 
layer (gray) and the Ti3d DOS of the interface Ti02 layer (orange), calculated both using 
LDA (Fig.[S4](a)) and LDA+£7 (Fig. [51(b)). Between 0.5eV and 2eV the main contributions 
to the total DOS of the full supercell arise from Ti3c?t 2 g orbitals located in the interface 
Ti02 layer. At energies below « 0.5 eV, 2p states located in the topmost AIO2 layer of 
the supercell contribute to the total DOS. 
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While Ti 3d states of the interface Ti02 layer almost completely dominate the calculated 
total DOS of the full supercell at energies between 0.5 eV and 2 eV and generate the three 
peaks between 0.6 eV and ~ leV, contributions from Ti3<i states of the adjacent Ti02 
layer considerably heighten the peak at 1.4 eV. In the total DOS of the full supercell (see 
Fig. IS4l) this peak appears to have a larger weight than found experimentally (Fig. 3(a)), 
which is possibly related to the current tunneling predominantly into electron states of the 
final Ti02 layer. 

Figure IS5l shows the orbital contributions to the Ti3<i DOS of the interface Ti02 layer 
both for LDA (Fig.ESKa)) and LDA+[7 (Fig.E3»). 

K. Calculated state densities of bulk LaAlC>3 and SrTiOs 

Figure EH shows the total state densities of bulk LaA10 3 (Fig.lSBTa)) and of bulk SrTi0 3 
(Fig.lSBTb)) calculated using LDA. The state densities are calculated employing the WIEN2k 
program package, using 3000 fc-points (84 fc-points in the irreducible wedge of the first 
Brillouin zone) for the LaAlOs DOS and 3375 fc-points (120 fc-points in the irreducible 
wedge of the first Brillouin zone) for the SrTiOs DOS. The DOS of the LaA103-SrTi03 
heterostructure is different from the DOS of either constituent. 
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FIG. SI: NDC(y s ) = (dI t /dV s )/(I t /V s + e)(Ks) characteristics with e = lpA/V acquired on sam- 
ples 1-4. The spectra were measured at 4.7 K with fixed tip-sample separations at sites far away 
from topographic steps. The different colors reflect different tip-sample separations. (a)&(b) Sta- 
bilization parameters: V^ab = 0.8 V, It,stab = 12 pA; lock-in parameters: V mo( i = 5mV rms , 
/mod = 77 Hz. (c) Stabilization parameters: Ki.stab = IV, it,stab = 12 pA; lock-in parameters: 
Vmod = 5mV rms , / mod = 77 Hz. (d) Stabilization parameters: T4 jSta b = 1 V, I tjS tab = 12 pA; lock-in 
parameters: V mo d = 10mV rms , / mo d = 161 Hz 6 . (e)&(f) Stabilization parameters: V^ab = 1.5 V, 
-^t.stab = 30 pA; lock-in parameters: V mo d = 2.83mV rms , / mo d = 977 Hz. (g)&(h) Stabilization 
parameters: V StSta b = 2.4 V, 7 t ,stab = 70 pA; lock-in parameters: V mod = 5mV rms , / mod = 77 Hz. 
The characteristics in (g) and (h) were measured at the same sample site. The spectrum shown 
in (a) is also plotted in Fig. 3(a) with a red line, the spectrum shown in (c) is also plotted in Fig. 
3(a) with a blue line, the spectrum shown in (h) is also plotted in Fig. 3(a) with a purple line. The 
data were averaged over an interval of 75 mV. As a guide to the eye, gray lines are drawn at the 
same positions as in Fig. 3(a). 
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FIG. S2: Comparison of the measured differential conductance and DFT calculated state densities, 
(a) NDC(V S ) = (dI t /dV s )/(It/V s + e)(V s ) characteristics with e = lpA/V measured at several 
sites located far away from topographic steps. The measurements were performed at 4.7 K with 
fixed tip-sample separations. The different colors reflect different tip-sample separations (blue: 
it,stab = 12 pA, y SiS tab = 1 V; red: / t ,stab = 12 pA, y SiSt ab = 0.8 V). The data were averaged over an 
interval of 75 mV. Both characteristics are also plotted in Fig. 3(a). (b) Ti3d DOS of the interface 
TiO"2 layer (orange) and O 2p DOS of the topmost AIO2 layer adjoining vacuum (green) calculated 
using LDA+t/ for a supercell with a 4 uc thick LaAlOs layer. In (a) and (b) the positions of 
characteristic features in the calculated Ti3<i DOS are marked with gray lines. 
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FIG. S3: Layer structure of the supercells used for the DFT calculations. 
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FIG. S4: Calculated total DOS of the full supercell with a 4 uc thick LaA10"3 layer on SrTiC>3 and 
state densities of the interface Ti02 layer. The total DOS of the full supercell (green), the DOS 
of the interface Ti02 layer (gray) and the Ti 3d DOS of this layer (orange) are shown. Note that 
the partial DOS (gray and orange) arise from integration over non-overlapping muffin-tin spheres, 
which neglects the remaining, interstitial regions. Therefore, the sum of the partial DOS is smaller 
than the total DOS (green), (a) State densities calculated using LDA. (b) State densities calculated 
using LDA+C7. 
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FIG. S5: Orbital contributions to DFT calculated Ti 3d state densities of the Ti ions in the interface 
TiO"2 layer. The total Ti 3d DOS (orange) contains contributions arising from Ti 3d xy orbitals (red) 
and Ti3d xz + 3d yz orbitals (blue), (a) Ti3d state densities calculated using LDA. (b) Ti3d state 
densities calculated using LDA+?7. 
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FIG. S6: State densities of bulk LaA103 (a) and bulk SrTi03 (b) calculated using LDA. 
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